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ABSTRACT 


This  research  contribution  presents  a  usage  manual  for  an  IBM  7090  computer 
program.  The  program  employs  a  Monte  Carlo  simulation  todetermine  the  prob¬ 
ability  of  destroying  individual  point  targets  within  a  target  complex  with  one  or 
more  groups  of  weapons.  It  is  assumed  that  the  groups  are  delivered  with  a 
bivariate  normal  aiming  error  and  that  the  individual  weapons  are  distributed 
with  an  independent  bivariate  normal  ballistic  dispersion.  The  program  is 
designed  for  conditional  damage  data  for  fragmentation  generated  by  an  IBM  7090 
program  furnished  by  the  U.S.  Naval  Ordnance  Test  Station  (NOTS),  China  Lake. 
A  flow  chart,  a  listing  of  the  FORTRAN  program  and  a  sample  problem  are  included. 
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I.  INTRODUCTION 


Previous  OKG  computer  programs  (reference  (a))  determine  the  probability 
of  destroying  rectangular  or  line  target s  with  a  single  group  of  weapons  delivered 
with  a  bivariate  normal  aiming  error,  anti  consider  the  individual  weapons  to  be 
distributed  about  their  mean  points  of  impact  within  the  group  with  a  bivariate 
normal  ballistic  dispersion.  This  program  determines  the  probabilities  of  des¬ 
troying  stationary  individual  point  targets  at  various  locations  within  a  target 
complex  with  one  or  more  groups  of  weapons  -  making  the  same  assumptions 
regarding  aiming  error  and  ballistic  dispersion  as  were  made  in  the  previous 
programs. 

In  considering  the  destruction  of  point  targets,  a  conditional  damage  function 
is  required.  The  conditional  damage  function  is  defined  in  reference  (b)  as  a 
function  which  gives  the  probability  that  a  target  which  is  located  at  a  particular 
range  and  bearing  from  the  weapon  detonation  point  will  suffer  at  least  the  stated 
degree  of  damage.  Analytical  functions  which  are  often  used  to  approximate  con¬ 
ditional  damage  functions  are  a  definite -range  (cookie -cutter)  function  or  a 
Gaussian  function  (these  functions  are  described  in  detail  in  reference  (b)).  Neither 
of  these  are  completely  adequate  however,  for  describing  the  effects  of  frag¬ 
mentation  on  a  point  target.  Therefore,  this  program  has  lieen  designed  to  use 
the  output  of  the  IBM  7U9U  warhead  lethality  program  written  at  the  U.S.  Naval 
Ordnance  Test  Station  (NOTS),  China  Lake.  The  N’OTS  program  provides, for 
each  5  degree  sector  around  the  weapon,  the  average  probabilities  of  destroying 
a  target  located  within  small  cells  at  various  radii  from  the  detonation  point.  It 
should  be  noted  that  the  design  of  the  program  does  not  preclude  the  use  of  a  cookie- 
cutter  or  Gaussian  conditional  damage  function  if  the  proliabilities  are  determined 
for  cells  at  the  proper  radii  (explained  in  a  later  section  of  this  research  con¬ 
tribution  for  each  5  degree  segment). 

Hi  GENERAL  DESCRIPTION 

Tins  program  is  a  Monte  Carlo  simulation  for  use  on  the  IBM  7090  computer, 
ft  determines  the  probabilities  of  destroying  individual  point  targets  within  a  target 
complex  on  the  basis  of  the  following  assumptions: 

a.  Tlie  aiming  error  and  the  ballistic  dispersion  are  bivariate  normal,  and" 
are  independent  in  the  same  coordinates,  hereafter  called  the  range  and  deflection 
coordinates  (the  names  applied  to  these  coordinates  result  from  a  common  appli¬ 
cation  of  the  model  to  bombing  problems  in  which  the  range  coordinate  is  in  the 
direction  of  the  flight  of  the  aircraft  and  the  deflection  coordinate  is  normal  to 
the  flight  path  in  the  horizontal  plane). 

b.  There  is  no  cumulative  damage  effect,  i.e.,  the  probability  of  destroying 
the  target  with  an  individual  weapon  is  independent  of  the  number  and  locations 

of  other  weapons  impacting  in  the  vicinity  of  the  target.  For  fragmentation  effects 
this  assumption  is  not  rigorously  valid  since  the  damage  criterion  is  often  stated 
as  some  minimum  number  of  fragments  per  unit  area  with  at  least  a  minimum 
energy  impacting  on  the  target.  Thus,  even  if  the  number  of  fragments  impacting 
on  the  target  from  any  one  of  the  weapons  might  be  less  than  the  number  specified 
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by  the  damage  criterion,  the  sum  of  the  fragments  from  two  or  more  weapons 
might  exceed  the  minimum  number. 

c.  The  conditional  damage  function  is  symmetrical  about  the  range  coordi¬ 
nate  through  the  weapon. 

Limiting  inputs  to  the  program  are: 

a.  Number  of  targets,  100. 

b.  Number  of  weapon  groups,  10. 

c.  Number  of  weapons  per  group,  100. 

d.  Number  of  conditional  damage  function  cells  per  5  degree  sector,  100. 

(Corresponding  to  a  maximum  radius  of  effect  of  3,920  feet). 

The  conditional  damage  data  must  lx-  provided  for  cells  at  various  radii  from 
the  weapon  for  each  5  degree  sector.  The  cells  -  except  for  the  innermost  -  are 
made  approximately  square  in  shape  by  selecting  the  outer  radius  of  each  cell 
to  be  equal  to  the  inner  radius  times  the  factor  (1  +  sin  5‘  ).  The  outer  radius  of 
the  inner  cell  for  each  sector  is  equal  to  1  foot  (see  figure  1).  The  value  of  the 
conditional  damage  function  is  taken  as  1  when  the  radius  is  less  than  1  foot, 
i.e.,  for  all  of  the  inner  cells.  The  individual  cells  are  denoted  by  a  pair  of 
subscripts,  k  I ,  where  k  refers  to  the  k^1  sector  numbering  from  the  tail  of  the 
weapon  and  1  refers  to  the  lta  cell  m  a  given  sector  numbering  out  from  the 
weapon.  The  program  accepts  the  NOTS  output,  in  which  angles  are  measured 
relative  to  the  bomb  tail  as  shown  in  figure  1 .  All  angles  calculated  in  the 
program  however,  are  measured  from  the  bomb  nose,  positive  direction  counter¬ 
clockwise;  but  because  of  symmetry  about  the  range  coordinate  of  the  bomb, 
values  of  the  conditional  damage  function  for  angles  that  are  measured  from  the 
bomb  tail  can  be  equated  by  a  simple  scheme  to  angles  that  are  measured  from 
the  bomb  nose;  i.e.,  a  cell  having  a  10  degree  angle  measured  from  the  bomb 
nose  would  have  the  same  conditional  damage  value  as  the  corresponding  cell 
in  figure  1  having  the  175  degree  angle. 

In  addition  to  the  parameters  specifying  the  weapon  effects  and  aiming  and 
ballistic  dispersions,  the  program  inputs  consist  of  parameters  specifying  the 
target  locations,  the  angle  of  approach  to  the  target  for  the  delivery  of  each 
group  of  weapons  (measured  from  the  X-axis  to  the  R-axis  -  see  figure  2),  the 
reliability  of  the  individual  weapons,  scale  factors  to  shrink  or  expand  the 
weapon  patterns,  the  number  of  dummy  passes  through  the  random  number 
generator  and  the  number  of  Monte  Carlo  iterations  to  be  made. 
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III.  MI-THUD  OF  SOLUTION 


Tin.  inputs  lux*  the  program  consist  of  tlu  following:  parameters: 


Add  re 

•ss’" 

Svmlxi! 

FORTRAN  Label 

Descript  ion 

1, 

1, 

1 

>'T 

ng 

NT 

Number  of  targets 

1, 

1, 

2 

NO 

Number  of  weapon  croups 

1, 

1. 

I 

i 

MCI 

Nunilx  r  of  Monti  Carlo  Herat  ions 

1, 

1. 

4 

n 

ni:p 

Numlx  r  of  empl\  passes  through 

r 

rjiulem  iiumlx  r  generator 

1. 

1. 

S 

Nc 

NC' 

Num.K  i*  of  damage  prohabilu\ 

com<  urs 

1, 

1, 

6 

.V 

ri;l 

Weapon  reliability 

1, 

1, 

7 

AR 

SlOR 

Range  aiming  error,  standard 
dev  tat  ion 

1, 

1. 

h 

AD 

SIGD 

Deflection  aiming  error,  standard 
de\ lat  ion 

1, 

!, 

4 

OR 

MR 

Range  ballistic  dispersion, 
standard  deviation 

1, 

1, 

10 

BD 

BD 

Deflection  ballistic  dispersion, 

standanl  deviation 

1, 

1, 

1  1 

gr 

GR 

Range  scale  factor  for  weapon 
group 

1, 

J  2 

gd 

GD 

Deflection  scale  factor  for 
weapon  group 

2. 

1...1, 

9 

j 

10 

N 

i 

NO) 

Number  of  weapons  in  i**1  group 

lf 

.1, 

I...1, 

L 

10 

t> . 

i 

M.PHA(I) 

Angle  of  approach  for  delivery 
of  ith  group  (degrees),  0  t>.<360 

1, 

4. 

1...  1, 

4, 

10 

AXi 

AX(  I) 

Aiming  point  of  i1^  group  X- 

coordinate 

1, 

5, 

1...1, 

5, 

10 

Av, 

AY(I) 

Aiming  point  oi  i1’1  group  Y- 

coordinate 

1, 

6. 

1 ...  1 , 

6, 

100 

X 

111 

X(M) 

Target  location  of  m1*1  target , 

X-cooi*dinate 

1, 

7, 

1 ...  1 , 7, 

too 

Yni 

Y<M) 

Target  location  of  m1"  target, 

Y-coordinat  e 

2, 

1, 

1...2, 

1,100 

ARij 

DELRO.J) 

Aimpoint  range  displacement  of 

m 

j1*1  weapon  of  i1*1  group  from  group 

• 

aiming  point 

4,  10. i... 2. 10, 100 


♦The  parameter  addresses  are  explained  in  section  FV,  User’s  Instructions. 


Address  Symbol  FORTRAN  Label  Description 

3. 1. 1. .  .3, 1, 100  Ap.  .  '  DELD(I,  J)  Aimpoint  deflection  displacement 

.  ^  of  weapon  of  i1*1  group  from 

;  group  aiming  point 

3.10.1.. . 3. 10. 100 

4. 1. 1 . .  .4, 1, 100  C, .  C(K,  L)  Probability  that  a  target  located 

.  in  the  k(h  sector  of  the  l1*1  annulus 

;  is  destroyed 


4,36, 1... 4, 36, 100 

Values  for  addresses  4, 1,  1  through  4, 36, 1  need  not  be  entered  by  the  user 
since  the  program  itself  sets  these  values  to  1.0.  The  user  may,  however,  enter 
whatever  values  he  wishes  for  these  cells  and  they  will  be  used  by  the  program. 

The  program  flow  chart  and  FORTRAN  statements  are  included  as  appendixes 
A  and  B.  The  solution  is  obtained  in  the  following  manner: 

1.  Before  starting  the  first  iteration  for  the  first  data  set,  n^  dummy  passes 

are  made  through  the  random  number  generator  (see  appendix  C).  Thus,  if  the 
same  data  are  run  more  than  one  time,  a  new  set  of  random  numbers  can  be 
selected  for  each  run. 


2.  For  each  weapon  in  each  group,  the  displacement  from  the  group 
reference  point  (A^,  Ap)  is  adjusted  with  the  scale  factors  (G^,  G^)  to  spread 

or  shrink  the  pattern: 


GR  V 


*D..  GDAD..  ^ 

ij  ij 

3.  For  each  Monte  Carlo  iteration: 

A.  Given  the  group  aimpoint  (A^.,  Ay.),  the  standard  deviations  of  aim 

error  (o^,  t™0  random  numbers  (n^,  n^)  selected  from  a  standard  nor¬ 

mal  distribution  (zero  mean  and  unit  variance)  and  the  approach  angle  o'.,  the 

coordinates  for  the  pattern  impact  reference  point  relative  to  the  X,  Y  axes  are 
determined  for  each  group: 

RXi=AXi  +  nl  CARcos‘'i+n2°ADsin‘'i  (3) 
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(10) 


and  for  0  s  R  <  1 ,  1  =  1 

R  2  1,  1  =  [4?.  --- :  —  +  2] 

ln(l  +  sin  5  )  I 

(where  [x]  is  defined  as  the  largest  integer  less  than  or  equal  to  x). 


F.  From  the  conditional  damage  data  for  the  kl^  cell,  gives  the 

probability  that  the  m^  target  is  destroyed  by  the  ij^  weapon.  Given  the 
reliability  (^),  the  probability  that  the  target  is  destroyed  by  at  least  one  of 
the  weapons  is: 


i  =  l 


O  -*ckl) 


(11) 


4.  The  Monte  Carlo  iteration  is  repeated  (i)  times  and  the  values  of  t 

are  accumulated  for  each  target.  At  the  conclusion,  the  estimate  of  the  proba¬ 
bility  of  destroying  each  target  is  determined  by: 


1 

T  i 

L*  m 

N=1 


(12) 


IV.  USER'S  INSTRUCTIONS 

An  OEG  Subroutine,  1-63S,  called  DATA  (see  appendix  E)  is  used  to  read 
the  punched  card  data  sets  into  the  computer.  A  feature  of  this  subroutine  is 
that  data  are  identified  and  stored  in  the  computer  memory  through  the  use  of 
addresses  that  are  the  subscripts  of  the  internal  data  array. 

The  addresses  shown  in  section  III  are  punched  into  cards  together  with 
the  values  of  the  parameters  stored  at  those  addresses.  A  data  form  which 
contains  the  addresses  and  parameter  values  is  shown  in  section  V  for  the  sample 
problem.  Any  number  of  successive  data  sets  may  be  submitted  at  one  time 
as  long  as  the  data  sets  are  separated  by  one  blank  card  and  the  last  data  set 
in  the  pack  is  followed  by  two  blank  cards.  Successive  data  sets  may  contain 
cards  for  only  those  parameters  that  have  values  different  from  the  preceding 
data  set. 


V.  SAMPLE  PROBLEM 

Four  targets  located  at  (1,  2),  (3,  -3),  (-1, 1)  and  (-2,  -1)  are  to  be  attacked 
with  two  groups  of  weapons  of  two  eacn.  Each  weapon  has  a  "cookie -cutter" 
damage  function  with  a  radius  of  effect  of  1 .  The  standard  deviations  of  aim 
error  are  3  in  tlv  range  coordinate  and  I  in  deflection.  The  ballistic  dispersion 
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(4) 


R,,.  =  A.,  -  n .  _ ,  n  sin  o .  -  n.,  ~ .  „  cos  o. 

'll  Yi  I  AR  1  2  AD  i 

(i  =1,2 . XG) 

11.  For  cadi  weapon  in  the  group,  the  group  impact  reference  point 
(R^.  Ry.),  the  weapon  displacement  from  the  reference  point  (j^.,  d|}1()- 

the  standard  deviations  of  ballistic  dispersion  (  "j^,  j^),  two  standard  normal 

random  numbers  (n^,  n^)  for  each  weapon,  and  the  approach  angle  e.,  the 

coordinates  of  the  weapon  impact  point  relative  to  the  X,  Y  axes  arc  determined: 


C.  For  each  target,  the  separation  distance  (squared)  from  each  weapon  is 
determined  using  the  target  location  (X  ,  Y  )  and  the  weapon  impact  points 

(X..,  Y  ): 
ij  U 

R2  =  (X  -  X  )2  +  (Y.  -  Y  )2  (7) 

ij  m  '  ij  m 

2 

D.  It  the  square  of  separation  distance,  R  ,  is  less  than  the  square  ol 
maximum  radius  of  effect  of  the  weapon,  R^  ,  the  relative  bearing  of  the 

target  from  the  weapon  is  determined  using  an  OFG  Subroutine,  ACOSD  (see 
appendix  D). 


=  l  cos 


.  /  X  -  X. 

•“  V" 


Y  -  Y.. 

m  ij 

'  MTV  -y.. 

\  u 


E.  The  conditional  damage  contour  cell  numbers,  k  and  1,  are  then 
determined  as  follows: 

for  0"  '' '  C  r  5!  ,  k  =  1 

5  -  1  s  10°,  k  =  2  (*>) 


175  M  x  180°,  k  =  36 
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i s  I  111  raiipi-  .mil  O.s  in  deflection.  TIk  ri  luhilitv  of  tin-  weapons  i"  O.s. 

I  In  centers  oi  both  groups  .in-  annul  at  tin  oriiini  ot  r hi-  \-\  a\i^  (H,d)  ami  t hi 
weapons  an-  spaced  alonp  tin-  raivac  coordinate  1  i  in  it  ^  apart.  One  itroup  is 
delivered  jt  an  aiiitlc  of  HI  to  tlii'  \-a.\is,  tin-  other  at  an  anjtle  of  I  Ml  to  t Ilf 
x-axis. 


FIG.  3:  SAMPLE  PROBLEM 

The  input  data  are  shown  on  the  attached  sample  data  form. 

VI.  KEYPUNCH  INSTRUCTIONS 

Keypunch  instructions  for  the  use  of  DATA  Subroutine,  I-63S,  are  con¬ 
tained  in  appendix  E. 

VII.  OPERATOR’S  INSTRUCTIONS 

Run  under  control  of  the  Bell  System  on  the  IBM  70d().  No  special  instructions. 


d 


VIII.  TIMING 


This  program  requires  approximately: 


1 

730 


x  N_  x  lx  r  N.  seconds 
1  i=l  1 


where 


N 

if 


N. 


l 


number  of  targets 

number  of  Monte  Carlo  iterations 

=  number  of  weapons  in  all  groups  summed. 
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i  .  ♦  r  l  a  G  :  I  INPUTS 


LI  .  ll  j 

tA*^  1  S 

«lAPCN  MlCl'P* 

iterations 

reliability 

random  numrer  masses 

, 

4 

2 

3000 

O.NCCO 

11* 

RANGE 

OEFLECTION 

AIMING  tKKOK  t  MEET 

3.0 

1.0 

BALLISTIC  iHRUR.  FEET 

1.0 

0.$ 

GROUP  SCALE  FACTOR 

1.0 

1.0 

GftCUP 

AIMING  POINTS 

APPROACH  ANGLE. 

NUMBER  OF 

NUKttER 

X  r 

DEGREES 

WEAPONS 

l 

1.  0. 

30.0 

2 

2 

C.  0. 

HO.O 

2 

GROUP 

WEAPCN 

AIMPCINT 

DISPLACEMENTS.  FEET 

NUMBER 

NUMBER 

RANGE 

DEFLECT  ION 

l 

I 

1.0 

0. 

l 

2 

-1.0 

0. 

2 

1 

1.0 

0- 

2 

2 

-1.0 
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CtbTKUCTION  PROBABILITY  CONTOURS 


K 

l 

1 

1 .  CGCCC 

2 

1  *  Llftlii! 

3 

1 . cccco 

4 

1  .  C  C  f  C  G 

5 

1 . CCCOO 

fc 

1 . COOCO 

f 

l . cccco 

a 

i .cccco 

y 

l. CCCOO 

10 

1 . COOOC 

li 

1 . occoo 

12 

1 . cccco 

13 

1 .cccco 

14 

1 . 00000 

15 

l. CCCOO 

16 

1 .cccco 

1  7 

1. CCCOO 

1  8 

l. 00000 

1  /■  • 

1  V 

l.OOCOO 

2  C 

1 .CCCOC 

21 

1. CCCOO 

22 

1.00000 

23 

l.OOCOO 

24 

l. cccco 

25 

l.COCOO 

26 

1 . 00000 

27 

l.OOCOO 

2b 

l.COCOO 

29 

1.00000 

30 

1.00000 

31 

l. occoo 

32 

l.COCOO 

33 

l.COCOO 

34 

l. 00000 

35 

l.COCOO 

36 

l.COCOO 

12 


M  ('I  -t 


TARGLT 

NUMBLR 


1 


X 

LCCAT lUNb 

y 

KILL 

PRLBAB I L I T  Y 

1  .0 

2  .0 

O.IIB 

3.0 

-3.C 

0 .  c  y 

-1.0 

1  .C 

0.^3^ 

-2.0 

-1 .0 

0.  184 
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CNA  COMPUTER  DATA  SUBMITTAL  FORM 

Submitted  by: 

Program  No. 


Special  Instructions: 


^/3-£>3  P  Fist.  Time  .  Classification  _ 


NOTES: 

1.  A  value  of  zero  must  be  entered  as  0,  not  left  blank. 

2.  Decimal  pts.  may  be  omitted  if  understood  to  follow  the  rightmost  digit. 

3.  The  value  3  X  10  ^  may  be  entered  as  .00003  or  3-5,  not  as  3  X  10 

4.  The  factor  portion  of  a  value  may  not  contain  more  than  8  digits. 

5.  The  exponent  portion  of  a  value  must  lie  within  the  range  ±39. 

6.  Exponents  may  be  omitted  if  zero.  If  not,  they  must  be  signed. 

7.  Blank  cards  should  be  indicated  by.  - - - b - . 
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References:  (a)  NAVWAG  Interim  Research  Memorandum,  NIRM-12;  "Usage 
Manual  for  a  Computer  Program  to  Compute  the  Effectiveness 
of  Groups  of  Weapons  against  Rectangular  and  Line  Targets" 
Unclassified  11  Dec  1959 

(b)  OEG  Study  626,  "Probability -of -Damage  Problems  of  Frequent 
Occurrence"  Unclassified  1 1  Dec  1959 
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FLOW  CHART 
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APPENDIX  B 


PROGRAM  LISTING 


I 


B-l 

(REVERSE  BLANK) 


DIMENSION  0(4,36.100).  DELR( 10, 100) .  DELO ( l 0 . 1 00  )  ,  P(IOO).  0(100). 
1T(100).  X(100).  Y(100).  N ( 10).  ALPHA (10).  AX(IO).  AY(10).  COSA(IO) 
2.  SINA(IO).  RN( 2 ) ,C  (  4  . 36 .100) . LL (10) 

EQUIVALENCE  (  O  .  C  )  • ( 0 ( 72 1 ) .REL  > • { D ( B65)  . S I GR  ) • (0( 1009 ) . S I  CO ) . 

1 (0(  1 153 ) .8H  )  .  ( D( 129  7 ) ,bO ) , ( D<  144 1  > . GR ) . ( D ( I 5AS ) . GO ) 

00  10  1=1.14400 
10  0(1 )=0.0 

DO  15  K= 1 »  36 
IS  C(4.K. 1 >=1.0 
REL  =  1,0 
GR=  1.0 
G0=  1.0 
NG=  1 

20  CALL  DATA  ( C  ,  4  ,  36 ,  1 0 0 .  ( NC ) 

IF  (INO-1)  30,420.20 

30  NT  =0 ( I  •  1  ,  1  ) 

NG=0( 1.1.2) 

MC 1 =0 (1.1,3) 

NEP  =  C(  1.1.4) 

NC=0( 1.1.5) 

ATEHM=1.0*SINDF<5.0> 

NCX  =  2  * ( NC- 1  ) 

HMAXSO  =  ATEI>M««NCX 
IF  ( NEP  )  41  . 4  1  ,  35 

35  CO  40  1=1. NEP 

40  CALL  GRNUMd  (CUMMY) 

011.1.4 ) =0. 0 

41  00  45  1=1. NG 

N<  I  )=D(  1,2.  I  ) 

ALPHA! I ) =01 1.3,1) 

AXl  I  J=0(  I  .4  ,  I  ) 

AY (  I  )  =  C(  1.5.  I  > 

I  I  =N  (  1  ) 

DO  45  J=1 . I  I 
OELRl l . J)=D<2, I.J) 

45  DELD1 I.J) =0(3, I.J) 

PRINT  50 

50  FORMAT  C 1 H I •  30X,  41HC  OMPLFX  TARGET  INPUTS) 

PRINT  51 

51  FORMAT  1////1JX,  8HC0N  TOURS «  5X.  7HTARGETS.  5X.  1 3HWEAP0N  GROUPS. 

1  5X ,  1  OH  I TERAT IONS.  5X,  1  1 HREL I AB I L  I  T Y ,  5X ,  20HRANDOM  NUMBER  PASS 

2ES  ) 

PRINT  52.  NC,  NT,  NG .  MCI.  REL.  NEP 

52  FORMAT  (//15X,  14.  9X ,  13,  12X.  12.  13X.  15,  IOX.  F6.4.  15X,  IS) 

PRINT  53 

53  FORMAT  ( ////42X.5HRANGE , 15X , 10HDEFLECT ION ) 

PRINT  54.  SIGR,  SIGO 

54  FORMAT  1//10X,  I  8HA l M I NG  ERROR.  FEF  T  ,  1 3X«  F6.1,  I7X.  F6.1) 

PRINT  55.  DR,  130 

55  FORMAT  (IOX,  2 1 HBALL I S T I C  ERROR.  FFFT.IOX.  F6.1.  1 7X.  F6.1) 

PRINT  56.  GR.  GD 

56  FORMAT  (IOX.  1 8HGR0UP  SCALE  FACTOR. 13X,  F6.1.  17 X,  F6.1) 

PRINT  57 

57  FORMAT  (//////12X,  5HGR0LP « 1 4  X ,  13HAIM1NG  POINTS, 1  IX.  15HA°PR0ACH 
1  ANGLE.. 9X,  9HNUMBER  OF) 

PRINT  50 

50  FORMAT  (  1  2X  ,  6HNUMUER.  11X,  1HX.  15X,  1HY,  13X,  7H0EGRFES.  I4X, 

1 7H*E  APONS ) 

PRINT  59 
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59  FORMAT  (lMO> 

DO  60  1=1. NG 

60  PRINT  61.  I.  AX||).  A  V  (  t  t  •  ALPHA ( I  ) .  Nil) 

61  FORMAT  (I3X.  12.  1IX.  F6.1.  tOX.  F6.I.  12X.  F5.I.  I6X.  13) 

PRINT  62 

62  FORMAT  I  t  H  |  .  12X.  SHGROUP  .  7X.  6HMt  APON  .  6X  .  28H  A  I  MPO I  NT  DISPLACEM 

IENTS.  FEET) 

PRINT  6  3 

6  3  FORMAT  (  I  3X  .  6MNUMHER.6X.6HNUMRER.9X , 5HRANGE . 7  X ,  I  OHOEFLF C T  I  ON  I 
PRINT  59 
00  6A  |=1 ,NG 
PRINT  59 
I  I  =  N  <  I  ) 

00  64  J= I  . I  I 

64  PRINT  65.1.  J.  DELHI!. J).  DFLO(I.J) 

65  FORMAT  (  1 4  X  *  12,  9  X ,  I  J.  9X,  F6.1,  9X.  F6.I) 

NCP -NC 

JL  =-  9 
JU  =  0 

67  IF  ( NCP- 10)  68. 68. 69 

68  JU=JU4NCP 
NCP  =  0 

GO  TC  70 

69  JU=JU4I0 
NCP-NCP- 1 0 

70  JL  =  JL ♦  I  0 
PRINT  71 

71  FORMAT  (1HI.  4  1  X  .  3? POL STRUCT  ION  t»l(OP  AH  1  L  I  I  V  CONTOURS) 

LL< I ) = JL 

CO  7?  1=2.10 

72  LLI  I  )  =  LL (  1-1)41 
LLL  =  JU- JL ♦ t 

PRINT  7  J,  I  LLI  I  ) ,  I  -  I  .LLL  I 

73  FORMAT  (  I  HQ  .  lux.  IHK,  f,X,  I0IIT.7X)) 

DO  74  K  K  =  1 ,  36 

74  PRINT  7  S  ,  KK.  I C I  4 . KK  . LNC )  .  LNC  =  JL.JU) 

75  FORMAT  I  1  OX  .  12.  6X,  I0IF7.S,  Jx)) 

IF  (NCP)  76,76.67 

76  N I  T  =  0 

CO  HO  I  =  I ,N  T 
80  Q(  I  >  =  0.0 

CO  90  1=1. NG 

COSA(l)=C OS CF( ALPHA!  1  )  ) 

SINAI  I  ) =SINCFI ALPHA l  I  )  ) 

JJ=N| I ) 

00  90  J=l  , JJ 
DELHI  I . J) =OELRI  I  ,  J ) *  OR 
90  OFLDI I ,J)=DFLOI I . J)«GC 
DO  95  M  = I .NT 
X  I M)=D(  I .6. M ) 

95  Y{ M)=D! I . 7, M) 

100  M  T  =  N 1 T ♦ 1 

DO  I  10  M= l , N  T 
I  l 0  T ( M )=  |  .0 
IG  =  0 

120  I G- I C  4 1 

DO  130  1=1.2 
130  CALL  CRNUMtJ  (  RN1  |  )  ) 
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RX  =  AX (  I  G  )  +RN  (  |  >»SIGR»C0SA<  I  G ) +RN( 21  *S  I GD*S 1 NA ( I G  ) 

RY  =  AV<  IG)*RN(|)»SIGH«SINA<  I  G 1-RN1 2) »S I GD»COSA ( IG> 

J  A  =  0 

140  J  V  —  JW ♦ 1 

IF  (BR+BD)  170. 170. ISO 
ISO  CO  160  1*1,2 
160  CALL  GRNUMB  (  RN (  1  )  ) 

170  XX=RX ♦ ( RN(  1  i •BR  +  DELR (  IG. JX  1  )  »C0SA(  I G ) ♦  (  HN ( 2 ) • BD* OELO (  IG. JN ) ) * 

1 SINAI IG) 

Y  Y  *RY  ♦  (  RN{  1  )  •0R4DELW  (  I  G  .  JH  )  )  *  S I  NA  (  I  G  >- (  RN  C  ?  )  •  BD4  DELD  <  IG.  J»  )  )  «COSA 
1  1  IG) 

M  T  =  0 

180  MT=MT*1 

IF  (T(MT))  300.300.iy0 
190  RSQ=(XX-X(MT))**2MYY-Y(MT))»*2 
IF  ( RSO-RMA X  SO )  200.200.300 
200  R=  SQRTF ( RSQ  ) 

IF  (R-I.O)  290.210.210 

2 1 0  THETA  =  ABSF(  ACOSDF  (  (X(Mf)  —  XX)/R)-SI GNF ( ALPHA (IG).Y(MT)-YY)) 

260  IF  < THETA-1  80.0 )  280,280.270 
270  THETa  =  360.0— THfcTA 

200  IF ( MODE! THcTA.  5.  ) )  202.281.282 

201  K=  THE  T  A/5 . 

GO  TO  283 

202  K=THETa/5.+  l  . 

203  L=LOGF(R1/LOGF( ATERM1E2.0 
K37=37-K 

T(MT)  =  T(MT)«(  1.0-C14.K37.L) *REL ) 

GO  TO  300 

290  T(MT)  =  T(MT)«(  1  .O-REl  ) 

300  IF  (MT-NT)  180,310,310 
310  IF  (JW-N(IGl)  140,320,320 
320  IF  (IG-NG)  120,330.330 
330  DO  340  M=  1 , NT 
340  0( M)*Qm+T(  M) 

IF  (NIT-MCI)  100,350.350 
350  DO  360  1 , NT 

360  P( M)*l  .0-01 M)/FL0ATF(MCI  ) 

PRINT  370 

370  FORMAT  1  1H I ,  12X,  6HTARGET.  I0X,  VHLOC A  T I ONS ,  14X,  4HKILL) 

PRINT  380 

300  FORMAT  (  1 3X ,  6HNUMBER ,  7X.  IHX,  13X,  1HY,  10X.  I  IMPROBABILITY) 

PRINT  59 
K=0 

DO  410  N  = 1 , NT 
K  =  K  ♦  1 

PRINT  390.  M,  X  C  M )  ,  Y  (  M )  »  P(M) 

390  FORMAT  (  1 3X ,  13,  6X,  FO.I.  0X,  F6.I,  7X,  F7.J) 

IF  { K— 1 0 )  410,400,400 
400  PRINT  59 
K  =  0 

410  CONTINUE 
GO  TO  20 
420  PRINT  430 
430  FORMAT  ( 1H1  ) 

CALL  ENCJOts 
END 
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APPENDIX  C 


GRNUMB  SUBROUTINE 


1.  Purpose: 

GRNUMB  provides  a  floating  point  pseudo-random  number  X  The  distribution 
of  successive  values  of  X  are  Gaussian  with  a  mean  of  zero  and  a  standard  devia¬ 
tion  of  one. 

2  Method: 

Consider  the  set  of  uniformly  distributed  pseudo -random  numbers  .  GRNUMB 
geneiates  a  sequence  of  Y.  by  the  method  of  congruences: 

Y.  =  2  35  (515235Yi_1  ,  mod  235) 

over  the  range  0  s  Y.  <  1  .  The  variance  of  this  uniform  set  is 

cY2  =  J  (Y- 1/2)2  dY  =  1/12  . 

o 

If  X  is  the  mean  of  any  selection  of  m  of  the  uniform  numbers  Y  ,  the  Central 
Limit  Theorem  states  that  the  variable  X  approaches  a  normal  distribution  where 
m  is  sufficiently  large.  A  satisfactory  value  for  m  is  30.  Values  of  X  are 
generated  as  a  sequence  of  Xn  ,  where  n  denotes  the  n^  entry  to  GRNUMB. 

. _  m  30 

Xn=Vl/m-Y2  l  (Y."1/2)  \  (Y.~1/2) 

i=l  i=l 

where  Y  =X  ,  ,  and  X  =  2'^.  The  variance  of  this  normal  set  is  1. 
o  n-1  o 

3.  Usage: 

X  is  obtained  by  use  of  the  statement: 

CALL  GRNUMB  (X) 

in  a  FORTRAN  program  for  the  IBM  7090. 

4.  Coding  Information: 

See  the  symbolic  listing  on  the  following  page.  GRNUMB  is  written  in  the 
7090  FAP  language  It  requires  40  words  storage  space  and  900  microseconds 
operating  time. 
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SYMBOLIC  LISTING 


GRNUMB 


C 


D 

E 


F 

XXI 

NUM 

NUMB 

MULT 

CHAR 

MAGIC 

H8 

H125 
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FAP 

REM  GRNUMB  G  WESTLUND  18  JUNE  1962  (7090> 

REM  GAUSSIAN  DISTRIBUTED  RANDOM  NUMBER  GENERATOR 

REM  ENTER  VIA  FORTRAN  STATEMENT  CALL  GRNUMB  (X: 

REM  SEQUENCE  STARTS  AT  DEC  1.  YIELDS  X  WITH  STAND.  DEV. 

ENTRY  GRNUMB 

SXA  XXL  l 

CLA  1.  4 

STA  F 

A  XT  30.  1 

STZ  NUM 

LDQ  NUMB 

MPY  MULT 

STQ  NUMB 

CLA  NUMB 

SUB  CHAR 

ARS  4 

ADD  NUM 

STO  NUM 

TIX  C,  1,  1 

LDQ  NUM 

MPY  MAGIC 

LRS  27 

TZE  D 

LRS  8 

CLA  H125 

ADD  H8 

LLS  8 

ALS  19 

TRA  E 

CLA  H125 

ALS  27 

STO  NUM 

CLA  H125 

LLS  27 

FAD  NUM 

STO  ** 

AXT  **,  1 
TRA  2.  4 
HTR  ** 

DEC  1 

DEC  30517578125 

TIX  0.  0.  0 

DEC  0. 31622780B0 

DEC  8 

DEC  125 

END 
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AS  IN  SUBROUTINE 

1 .  Purpose: 

Given  the  floating  point  argument  X,  ASIN  provides  a  floating  point  number 
Y,  where 

Y  =  sin  ^X  or  Y  =  cos^X 

-I  ?  X  r  1  -U-Xsl 

-  ”/2  SYS  -/2  OtYs" 


Depending  upon  the  entry  to  the  subroutine,  Y  may  be  found  in  radians 
or  degrees . 


2.  Method: 

A  series  approximation  (Hastings),  is  used  to  find  the  Ar~  cosine.  If  the 
Arc  sine  is  desired,  the  arc  cosine  is  subtracted  from  r/ 2.  The  result  is  then 
converted  from  radians  to  degrees  if  desired. 


Arcsin  X  =  ^  -  /1-X  ^(X) 

$(X)  =  ao  +  a^X  +  a2X  +  •  •  ♦  +  a  yX 


a0  =  1.5707,  9630,  50 
=  .2145,  9880,  16 
a2  =  .0889,  7898,  74 
a3  =  -.0501,  7430,  46 


a4  =  .0308,  9188,  10 

a.  =*.0170,  8812,  56 

a,  =  .0066,  7009,  01 
o 

ay=  -.0012,  6249,  11 


degrees  =  57.2957,  7951,  3  radians 


3.  Usage: 

Y  is  obtained  by  use  of  one  of  the  four  expressions  below  in  a  FORTRAN 
program  for  the  IBM  7090. 


ASINF(X) 

ACOSF(X) 

ASENDF(X) 

ACOSDF(X) 


Arc  sine  X  (radians) 
Arc  cosine  X  (radians) 
Arc  sine  X  (degrees) 
Arc  cosine  X  (degrees) 


4.  Coding  Information: 

See  the  symbolic  listing  on  page  D-2.  ASIN  is  written  in  the  7090  FAP 
language.  It  requires  51  words  storage  space  and  850  microseconds  operating  time. 
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SYMBOLIC  LISTING 


ASIND 

ACOSD 

AS  IN 
ACOS 

A 

B 

C 

D 

E 

X4 

S 

X 


FAP 

REM  3-62S  ASIN  G.  WESTLUND  7090 

REM  ARCSINE -ARCCOSINE  SUBROUTINE.  ENTRY  VIA  FORTRAN 
REM  EXPRESSION  ASIN(X)  OR  ACOS(X)  YIELDS  RADIANS. 

REM  ASIND(X)  OR  ACOSD(X)  YIELDS  DEGREES. 

REM  ARCSINE  RANGE  =  -90  to +90,  ARCCOSINE  RANGE  =  0  TO  180. 

ENTRY  ASIN,  ACOS,  ASIND,  ACOSD 

SXD  S,  4 

SXA  S,  4 

TRA  *+2 

SXD  S,  4 

SXA  X4,  4 

STO  SIGN 

SSP 

TZE  D 

STO  X 

CLA  ONE 

FSB  X 

TZE  e 

CALL  SQRT 

STO  R 

AXT  7,  4 

LDQ  T+7- 

FMF  X 

FADT+7,  4 

XCA 

TIX  B,  4,  1 
FMP  R 
FAD  PI  2 
LDQ  SIGN 
LLS  0 
LXD  S,  4 
TXH  E,  4,  0 
CHS 
FAD  PI2 
LXA  S,  4 
TXL  X4,  4,  0 
XCA 

FMP  CONV 
AXT  **,  4 
STZ  S 
TRA  1,  4 
HTR  O 
HTR  ** 
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aS  H 


SYMBOLIC  LISTING  (Cont'd) 


SIGN  HTR  ** 

ONE  DEC  1 . 

PI2  DEC  1 . 5707963 

HTR  ** 

DEC  -.0066700901,  .0170881256,  -.0308918810,  .0501743046 
DEC -.0889789874,  .2145988016,  -1.570796305,  .0012624911 
CONV  DEC  57.295779513 

END 
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APPENDIX  E 


DATA  SUBROUTINE 


!  .  Introduction: 

Many  computer  programs  require  the  flexibility  of  varying  any  or  all  of  the 
parameters  in  a  computer  run.  Although  FORTRAN  is  fairly  flexible  in  its 
arithmetic  and  control  statements,  its  input-output  statements  are  .juite  rigid. 

In  order  to  read  cards  for  instance,  considerable  effort  must  be  expended  by  the 
FORTRAN  programmer  in  writing  his  input  statements.  This  subroutine  eliminates 
some  of  that  tedium.  The  concept  of  a  "data  set"  is  used.  A  data  set  consists  of  a 
sequence  of  punched  cards  terminated  by  one  blank  card.  A  parameter  deck  for  a 
computer  run  may  consist  of  several  data  sets.  Such  a  parameter  deck  is  terminated 
by  two  blank  cards. 

2.  Parameter  Addresses: 

The  primary  advantage  of  this  subroutine  over  FORTRAN  input  statements  re¬ 
sults  from  the  use  of  "parameter  addresses."  An  address  is  a  relative  location  in 
the  computer  memory.  It  is  the  subscript  of  an  array  or  matrix.  For  example,  in 
an  array  called  X,  the  parameter  value  would  be  located  at  address  53.  By 

using  the  parameter  addresses,  a  user  of  the  program  need  submit  only  those  param¬ 
eter  values  in  a  data  set  that  are  different  from  those  in  the  previous  set. 

Three  types  of  addresses  are  permitted  by  this  subroutine. 

(1)  A  numeric  address  consisting  of  one  to  five  characters,  each  of  which 
is  a  digit  0-9.  Such  an  address  (n)  refers  to  the  nth  element  in  a 
specified  array. 

(2)  An  alpha  address  consisting  of  one  to  six  characters,  the  first  of  which 
must  be  alphabetic  (A-Z).  The  remaining  may  be  alphabetic  or  numeric 
(A-Z  or  0-9).  Such  an  address  refers  to  the  n”1  element  in  a  specified 
array  (1  s  n  <  26),  where  the  first  character  of  the  address  corresponds 
to  n  as  the  26  letters  of  the  alphabet  correspond  to  the  integers  1-26. 

(3)  A  matrix  address  consisting  of  two  or  more  numeric  fields  separated  by 
commas.  For  example,  the  address  53,  47  refers  to  the  element  in  the 
53"  *  row  and  the  47t'>  column  of  a  two-dimensional  matrix.  There  is  no 
Unui  to  the  number  of  dimensions  in  a  matrix  address. 

3.  Input  Card  Format: 

A  standard  submittal  form  (see  attachment)  has  been  designed  for  the  analyst. 
This  form  provides  for  entering  parameter  values  with  their  associated  addresses. 
The  user  indicates  blank  cards  to  separate  data  sets.  The  keypunch  operator  has 
the  option  of  punching  one  address  and  value  per  card,  or,  if  the  addresses  are 
sequential,  of  punching  one  address  and  several  values  on  a  card. 
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Only  columns  1-72  of  a  card  are  used.  Each  column  nrv  st  contain  one  of  the 
following:  a  digit  (0.-9),  a  "  +  "  or  "  -  "  sign  or  a  dash,  a  letter  (A-Z),  a  period, 
a  comma,  or  a  blank.  Each  punched  card  must  contain  one  parameter  address. 

The  address  may  start  in  column  1.  or,  if  des«red,  may  start  in  a  later  column, 
provided  all  columns  before  it  are  blank.  The  address  is  terminated  by  at  least 
one  blank  column.  Only  one  addresa  is  permitted  on  the  card.  Succeeding  col¬ 
umns  contain  one  or  more  parameter  values,  each  separated  by  one  or  more 
blank  columns,  A  value  may  be  signed  or  unsigned.  The  length  of  the  value  field 
in  variable.  No  blanks  are  permitted  within  a  value  field.  A  value  may  be  punched 
with  or  without  an  exponent.  An  exponent  is  recognized  by  the  presence  of  a  plus 
or  minus  sign  (or  dash)  between  the  fractional  part  and  exponent  part  of  the  value. 
Decimal  points  (periods)  may  be  punched  in  either  the  fractional  or  exponent  parts 
of  a  value.  If  more  than  one  value  is  punched  on  a  card,  those  after  the  first  will 
be  entered  at  sequential  addresses  relative  to  the  address  of  the  first  value. 

4.  Usage: 

A  data  set  is  read  by  the  use  of  the  statement: 

CALL  DATA  (X,  I) 

in  a  FORTRAN  program  for  the  IBM  7090.  The  argument  X  is  the  name  of  an 
array  in  the  program.  The  argument  I  is  an  indicator  set  by  the  subroutine. 

This  indicator  may  be  tested  by  the  main  program  upon  return  from  the  sub¬ 
routine.  It  will  have  a  value  of  0  or  1  or  2  . 

0:  The  subroutine  has  read  a  data  set.  The  main  program  will  normally 

proceed  to  operate  on  this  data. 

1:  The  subroutine  has  read  the  second  blank  card  which  terminates  the 

parameter  deck.  The  main  program  will  normally  terminate  at  this 
point. 

2:  The  subroutine  has  read  a  "bad"  data  card.  The  main  program  may 

terminate  the  run,  or  ignore  the  card  and  return  to  the  subroutine  to 
read  the  rest  of  the  data  set. 

If  the  cards  to  be  read  contain  matrix  addresses,  additional  arguments  must  be 
included  in  the  FORTRAN  calling  statement: 

CALL  DATA  (X,  D. ,  D0.  D~,  .  . . ,  D  ,  I) 
i  i  o  n 

where  D.  is  the  i1^  dimension  of  the  matrix  X. 

5.  Method: 

See  the  attached  flow  chart.  DATA  reads  parameter  values  and  loading  addres¬ 
ses  from  cards.  If  sense  switch  5  is  up,  it  will  read  the  values  and  addresses  from 
tape  (unit  A2).  It  converts  the  values  to  floating  point  numbers,  and  stores  them  as 
elements  of  an  array  specified  in  the  calling  statement.  The  elements  are  specified 
by  the  addresses.  If  a  card  (or  tape  record)  is  read  which  contains  non-permitted 
characters  (see  input  card  format  above),  DATA  prints  the  statement  "bad  data  card," 
followed  by  an  image  of  the  card  itself. 
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